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We present a semi-numerical method to compute one-loop corrections to processes involving many 
particles. We treat in detail cases with up to five external legs and massless internal propagators, 
although the method is more general. 
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... I. INTRODUCTION 
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' The importance of next-to-leading order (NLO) corrections is today well-established. The benefits of NLO compu- 
tations include the possibility of reliably estimating cross-section normalizations, of reducing renormalization scale 
K*" . dependences and of understanding the uncertainties due to the perturbative expansions. For all searches of new 
physics at upcoming collider experiments it is crucial to know the background in good detail, the more so, the 
smaller the signal to background ratio is. At the LHC and at the ILC most processes will involve many particles in 
the final state. It is therefore of increasing importance to predict cross sections for those processes at NLO. In QCD 
all 2 ^ 2 processes are today known, however a few 2-^3, most 2^4 and all 2 ^ processes with N > 5 are not 
yet known at NLO. 

A full A^ particle NLO calculation requires i) a tree level (N-l-l)-particle amplitude; ii) the NLO N-particle ampu- 
le , tude; iii) the subtraction terms to regulate the divergences of both, i) and ii). While point i) has been extensively 
CO I treated in the literature and automatized 0], and point iii) is also well- understood the bottleneck in N-particle 
I ' NLO calculations remains the complexity of the analytical evaluation of the virtual contribution. The aim of this 
i project is therefore to seek a semi-numerical solution to this. Here we report on early steps along this direction. For 
other recent progress in developing algorithms to evaluate one- loop integrals see for instance refs. 

^ ■ 

p • II. SEMI-NUMERICAL METHOD 

(D : 

'T^ We first introduce some notation. We define the D dimensional N-particle M-tensor integral as 

^> : 

k>( ' f d^l l^^ l^^" 

^ , J ITT / (Ji . . . fZjY 

Notice that we consider here only the case with massless internal propagators, although the method is more general. 

A. The algorithm 

The numerical procedure we used is simple and is based on the following few steps 

1. use a Feynman diagram generator (e. g. Qgraf |^) to generate the amplitude A for a specific process at NLO; 

2. use a symbolic manipulation program (e. g. Form 3|) to write the amplitude as 

A{pi,...,pN;ei, ...eN;...) = ^K^^,,,^^{pi, . . . ,pn; ei, . . . eN; ■ ■ .) ■ I^^'-'^" {D;vi, ...,vn), (2) 

n 

where the kinematic tensor K depends on the particle properties (momenta pi, polarization e^, . . . ) and is 
made up of four-dimensional objects only (i. e. dependence on the metric musthe canceled analytically); 
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3. use Davydychev reduction formula 

/Mi...MMp;W---,^iv})- E (-1) M^i^i]^' ■■■[qNm^,...^,A'^i)n.---M.. (3) 



X /(D + 2(M- A);j/i +Ki,...,i/Ar + K7v), 



to reduce any tensor integral to higher dimensional scalar integrals; 
4. use the basic equation of the integration-by-parts method 



d 



= 0, (4) 



N 



to derive a complete set of reduction relations. Here complete means that any integral is reduced to a linear 
combination of analytically known scalar integrals. A sample reduction relation one obtains is for instance 



(5) 



where = {qi - q^Y , = Sij ^ and B = YliLi 



The method is semi- numerical in the sense that point 1) and 2) are done analytically, once and for all for a given 
process, while step 3) and 4) are repeated numerically for each phase space point. A key point which makes this 
method efficient is that a record is kept of all previously computed scalar integrals, so that each one is computed 
only once. 

B. Treatment of exceptional phase space points 

Implicit in the numerical use of recursive relations such as eq. © is the assumption that the kinematic matrix 
Sij is not singular and that B does not vanish. If B or det(S') are exactly zero one obtains a simpler set of 
relations Q. This is needed for instance in the treatment of cases with high N, where the particle momenta are not 
all linearly independent. More problematic is the treatment of so called exceptional momentum configurations, where 
an accidental degeneracy causes B or det(S') to be very small. In this case, relations such as ec^. ^ are in principle 
still valid but become numerically unstable. A solution to this problem was first suggested in [lOj. The idea is to 
exploit the existence of a small quantity, parameterizing the closeness to the exceptional phase space configuration, 
to define expanded reduction relations. For instance if B <C 1 one rewrites eq. (0 as 

N 

m = E ^'^(^; {^fc - + {D + l-a)BI{D + 2; {^.4^1) , (6) 

1=1 

where the first term consists of simpler integrals (lower D and/or a) and the second term consists of more difficult 
integrals, which are however suppressed by the small parameter B in front. In a similar way, if det(S') is small one 
determines the eigenvalue(s) corresponding to a small eigenvector and defines modified, expanded reduction relations. 



C. Higgs plus four partons 

As a first application of the method we considered the yet uncalculated gluon-gluon fusion amplitude pi + p-^ —> 
Pd.+Pi + H aX NLO in the large mt limit, where the Higgs couples directly to the gluons via an effective coupling 1 1] . 
The main mechanisms for Higgs production at hadron colliders are vector boson fusion (VBF) and gluon gluon fusion 
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(GGF). The former allows a precise determination of the Higgs properties, in particular of the Higgs couplings. Due 
to the gluons in the t-channel, GGF is characterized by more QCD activity in the central region, therefore VBF and 
GGF can be partly discriminated with suitable kinematical cuts. However, GGF remains the dominant background 
to VBF. A precise determination of it is therefore mandatory. We performed the numerical, partonic calculation E3 |. 
In the four quark case we performed also an analytical calculation |l2j , which confirmed the validity of the numerical 
results. In the two quark-two gluon and in the four gluon case we verified the Ward identities and known identities 
between gluonic amplitudes (cyclicity, reflection and dual Ward identity and the decoupling identity for number of 
flavours Uf = 0) ^] . Since the method is based on an analytical evaluation of basis integrals no loss of accuracy is 
expected. Using double precision the accuracy of the results is of the order of 10^^^ (10^^) for normal (exceptional) 
phase space points (a even higher accuracy can be achieved at the price of including more terms in the expansions) . 

III. WORK IN PROGRESS 

The combination of real and virtual results and phase space integration to obtain full predictions for Higgs + dijet 
production at the LHC is currently begin performed. We are also exploring the possibility of refining the method 
and of finding a more efficient numerical algorithm and phase space integration procedure. 
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